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The LHCb collaboration has recently reported on some anomalies in b —► s transitions. In addition to dis¬ 
crepancies with the Standard Model (SM) predictions in some angular observables and branching ratios, an 
intriguing hint for lepton universality violation was found. Here we propose a simple model that extends the 
SM with a dark sector charged under an additional t/(l) gauge symmetry. The spontaneous breaking of this 
symmetry gives rise to a massive 7! boson, which communicates the SM particles with a valid dark matter 
candidate, while solving the > s anomalies with contributions to the relevant observables. 


I. INTRODUCTION 


Rare decays stand among the most powerful probes of 
physics beyond the Standard Model (SM). Indeed, most new 
physics (NP) scenarios suffer from severe constraints due to 
their potentially large contributions to flavor observables. This 
has motivated an intense experimental search, with special fo¬ 
cus on observables which are strongly suppressed in the SM. 

In 2013, the LHCb collaboration reported on the measure¬ 
ment of several observables in processes involving b—>s tran¬ 
sitions. A significant tension with the SM was found in some 
cases. These include angular observables in B —)■ K /j p (T|, 
particularly large in case of the popular p\ GUI, as well as 
a decrease, with respect to the SM expectation, in several 
branching ratios j5j |6). These anomalies received immedi¬ 
ate attention in the flavor community, and soon several inde¬ 
pendent global fits QUID showed that the tension could be 
alleviated in the presence of new physics contributions. Inter¬ 
estingly, in 2014 the LHCb collaboration also found an indica¬ 
tion of lepton universality violation in the theoretically rather 
clean ratio hd 


BR (B^Kfi+p-) 

BR (B-A-Ke+e-) 


0.745 


+ 0.090 

- 0.074 


±0.036, 


( 1 ) 


measured in the low dilepton invariant mass regime, which de¬ 
parts from the SM result R| m = 1.0003 ±0.0001 by 2.6a lfL2l 
Again, this led to some excitement in the community, in par¬ 
ticular after it was found that this hint is compatible with the 
previous anomalies in b —> s transitions: they can be explained 
by the same type of NP contributions to muonic operators C3- 
m So far, the discussion has focused on results obtained 
by LHCb with an integrated luminosity of 1 tb 1 . The lat¬ 
est chapter of this story is the recent announcement of the 
LHCb collaboration of new results using the full LHC Run 
I dataset lfl8ll . with an integrated luminosity of 3 fb 1 , which 
confirm the robustness of the LHCb data. Indeed, the new re¬ 
sults are compatible with those found with 1 tb , and several 
theorists have used the new data to update their analyses. Ac¬ 
cording to HD, the hypothesis of non-zero NP contributions 
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is preferred over the SM by 3.7c, 4.3c> if the 2014 measure¬ 
ment of Rk is included, whereas the analysis of l20l finds a 
slightly larger statistical significance for NP (slightly larger 
than 4a) even in the absence of Rk- In any case, assuming 
that hadronic effects EDE) are not behind these anomalies 
(something impossible in case of Rk), there is clear evidence 
of new physics in B meson decays. For a complete review of 
the subject see 02) and references therein. 

Several approaches have been considered in order to ex¬ 
plain the b —> s anomalies. Some papers 12314251 consider 
generic Z' bosons with flavor violating couplings. These have 
been shown to provide a simple way to reconcile theory pre¬ 
dictions with experimental data. There are also some works 
that were built to address the first anomalies but fail to ad¬ 
dress the lepton universality violating Rk measurement, see 
for example Refs. |26l|27) in the context of 331 scenarios 
or Ref. fl28l . where 4-quark scalar interactions were consid¬ 
ered. Finally, a few recent models can also account for the 
lepton universality violating Rk measurement. One finds three 
types of working models: models with scalar or vector lepto- 
quarks (or similarly, R-parity violating supersymmetry, where 
the squarks can play the role of the leptoquarks) 1141129143 11 . 
composite Higgs models l25][32l and models with a Z!. The 
latter is generally considered the easiest approach, since the 
anomalies can be solved with new contributions to vectorial 
operators. However, to the best of our knowledge, the only 
(complete) models that have been put forward in order to ac¬ 
count for these anomalies using a Z' are: (1) the one intro¬ 
duced in Ref. [f33l (and the Two-Higgs-Doublet version in 
<341 ). which makes use of a U (1 )/ Jfl -+ T gauge symmetry, and 
(2) the one proposed in the recent paper f35l . which extends 
the horizontal gauge symmetry to the quark sector. 

The existence of dark matter (DM) is a well established evi¬ 
dence of new physics, supported by a plethora of astrophysical 
and cosmological observations. This motivates further exten¬ 
sions of the SM of particle physics that include valid DM can¬ 
didates. In DM models, the current Planck 3a limits for the 
DM relic density, 0.1118 < mIi 2 < 0.1280 lf36l . typically 
imply strong constraints on the mass and couplings of the DM 
candidate. This, combined with additional constraints (such 
as those from flavor physics), sometimes leads to very pre¬ 
dictive scenarios where the parameter space of the DM model 
shrinks to small regions compatible with all observations. 

We propose a simple model which captures the main in- 




2 


gredients required to explain the b —7 s anomalies and, simul¬ 
taneously, provides a DM candidate. In order to do so, we 
extend the SM with a dark sector charged under an additional 
1/(1) gauge symmetry. The spontaneous breaking of this sym¬ 
metry gives rise to a massive 2 ! boson, which communicates 
the SM particles with the dark matter particle and solves the 
b —7 s anomalies with contributions to the relevant flavor ob¬ 
servables. The interplay between DM and flavor physics leads 
to a very constrained scenario. In particular, achieving the 
measured DM relic density while providing the required NP 
contributions to B meson decays turns out to be very restric¬ 
tive on the model parameters. 

The rest of the paper is organized as follows: in Sec. [H|we 
discuss some general aspects about the b —7 s anomalies and 
introduce the basic language to be used throughout the paper. 
In Sec. m we define our model, show how it addresses the 
h -7 s anomalies while providing a DM candidate and discuss 
extensions to generate non-zero neutrino masses. Sec. TV| is 


devoted to the DM phenomenology of our model and the Z' 
portal responsible for its production in the early universe. Sec. 
[V|is a review of the most relevant constraints in our model, to 
be considered in Sec. m where we present our numerical 
results. Finally, in Sec. VII we summarize our results, discuss 
some related aspects and derive some general conclusions. 


to the unprimed Wilson coefficients, which can be written as 
C 9 = C 9 SM + Cf\ (5) 

C 10 = Cf 0 M + Cr 0 P . (6) 

The SM contributions, Cf M and C™, have been computed by 
different groups. Assuming that these are the only contribu¬ 
tions to the Wilson coefficients, several independent global fits 
have found a sizable tension with experimental data on b — 7 s 
transitions. This motivates the addition of NP contributions. 
When this is done, the global fits are clearly improved. Ac¬ 
cording to ma (we will mainly consider the results of this 
global fit), the best improvements are found in two cases: 

• Scenario 1: NP provides a negative contribution to 0 9 , 
with Cg' NP ~ —30% x cf M , leading to a Wilson coef¬ 
ficient Cg significantly smaller than the one in the SM. 

• Scenario 2: NP enters in the SU (2 ) R invariant direction 
C*’ NP = -C^ NP , with Cf v ~ -12% x cf M . 

In both cases, the rest of operators involving muons are per¬ 
fectly compatible with the SM expectations. Similarly, no NP 
is required for operators involving electrons or tau leptons. 


II. THE b > s ANOMALIES 

In this Section we discuss some general aspects about the 
b -7 s anomalies. For a complete review of the subject, we 
refer the reader to ini. 


B. Model building requirements 

Once determined the type of contributions to Cg and Cjo a 
NP model has to induce, one can figure out the main ingredi¬ 
ents of a simple working model. Arguably, the simplest one 
contains the following elements: 


A. Operators and global fits 


• A massive Z 1 boson, responsible for the vectorial oper¬ 
ators Og and Oio 


The effective Hamiltonian for b — 7 s transitions is usually 
written as 

4Gf e 2 

V,bV* ^ £ (QOi + c; O') + h.c., (2) 


• The Z' must have flavor violating couplings to quarks 

• The Z' must couple differently to electrons and muons 
This setup can be easily parameterized by the Lagrangian 

HUE) 


where Gp is the Fermi constant, e the electric charge and V 
the CKM matrix. 0, and O' are the effective operators that 
contribute to b — 7 s transitions, and C, and C\ their Wilson 
coefficients. Since the most important anomalies have been 
found in semileptonic B meson decays, we will consider the 
following set of operators, 

Og = (sy,P L b) {efi) , O’g = {s^P R b) {ere) , (3) 

Oio = {*r p L b ) {bflsi) , 0 [o = {sy„P R b) {If^e) . (4) 

Here £ = xQ It is also customary to split the Wilson co¬ 
efficients in two parts: the SM contributions and the NP con¬ 
tributions. Since the primed operators, 0' 9 and O[ 0 , do not re¬ 
ceive significant SM contributions, this is usually applied only 


L D 


fit (A 1 U Pl 


+ A f / i P R )f j Z' tI . 


(7) 


and either ( 1 ) A^“ = 7 ^ 0 , or ( 2 ) 7 ^ 0 and A 1 )) 1 = 0 , 


In order to account for the anomalies, one requires A( M 7 ^ 0 


depending on the scenario one wants to consider. The rest of 
the Z' couplings to SM fermions can be set to zero. As we will 


see in Sec. Ill the model we are going to consider belongs 
to scenario (2). In this case, the quark and lepton currents 
contributing to Og and Oio are both left-handed and one finds 
at tree-level fT7l 


rftNP _ ^,NP _ A, A l 

~ L 10 — v v * 

Tf/j Us 


m Z ' 


( 8 ) 


with 


1 When referring to operators involving a particular lepton flavor, we will 
denote it with a superscript, e.g. Cg and Og, for muons. 


Av 


, c y/ 2 

\ y/lGpOL) 
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4.94 TeV, 


(9) 
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where a = £-7 is the fine structure constant. Note that A v and 
the CKM elements appear in Eq. ([8]> in order to normalize the 
Wilson coefficients as defined in Eqs. (|3]i and (|4j. 

We now introduce a complete renormalizable model with 
these properties. 

III. THE MODEL 

We extend the SM gauge group with a new dark U(\)x fac¬ 
tor, under which all the SM particles are assumed to be sin¬ 
glets. The only particles charged under the U(l)x group are 
the following vector-like fermions, 

Ql = ^3,2, ^,2^ , 3,2, i,2j , (10) 

L l = , L R =(l,2,~,2) , (11) 

as well as the complex scalar fields 

4 >= ( 1 , 1 ,0,2), x = ( 1 , 1 ,0,-1), (12) 

where we denote the gauge charges under SU (3) c ® SU(2)l ® 
U(l)y <8> U(l)x and the SU(2) R doublets can be decomposed 
as Ql,r = ( U,D) lr and L l r = (N,E) LR . 

Besides canonical kinetic terms, the new vector-like 
fermions have Dirac mass terms, 

L m =m Q QQ + m L LL, (13) 

as well as Yukawa couplings with the SM fermions 

L Y = XQQ R <\>q L + \ L L R (\)£ L + h.c ., (14) 

where Xq and X/ are 3 component vectors. The scalar poten¬ 
tial takes the form 

v=v SM + r(H,^x) + r(^x). (is) 

Here H is the SM Higgs doublet and Tsm = + \\H \ 4 

is the SM scalar potential. The pieces involving the U(l)x 
charged scalars are 

= X H * \H \ 2 \$\ 2 +X Hx \H\ 2 \x\ 2 (16) 

and 

^(<t),x)=^l<l>| 2 +^|x| 2 + ^|(|)| 4 + ^|x| 4 

+^zl^| 2 |x| 2 +(^te 2 + h.c.) . (17) 

We will assume that the scalar potential is such that only the 
standard Higgs boson and the <f> field acquire non-zero vacuum 
expectation values (VEVs), 

(H°) = ^=, <*|>) = (18) 


where gx is the U (\)x gauge coupling, and inducing mix¬ 
ings between the vector-like fermions and their SM counter¬ 
parts thanks to the Yukawa interactions in Eq. ©■ Further¬ 
more, after spontaneous symmetry breaking, the resulting La- 
grangian contains a remnant Z 2 symmetry, under which x is 
odd and all the other fields are even. Therefore, X is a stable 
neutral scalar, and thus a potentially valid DM candidate. It is 
worth noting that the mechanism to stabilize the DM particle 
does not introduce additional ad-hoc symmetries, but simply 
makes use of the same U(l)x symmetry that is required in or¬ 
der to give an explanation to the LHCb observations. This 
goal has been achieved by breaking the continuous U(l)x 
symmetry to a remnant Z 2 , something that can be easily ac¬ 
complished with a proper choice of U (\)x charges 13814401 . 

Before concluding this section we must comment on (7(1) 
mixing. It is well known that nothing prevents 1/(1) factors 
from mixing. In the model under consideration, this would be 
given by the Lagrangian term ATI 

LDeF*F?, (19) 

where K are the usual field strength tensors for the U(\)x,y 
groups. In the presence of a non-zero £ coupling, kinetic mix¬ 
ing between the t/(l)x and U(l)y gauge bosons is induced. 
As a consequence of this, the physical Z! boson would cou¬ 
ple to all particles that carry hypercharge, this is, to all the 
SM fermions. This would lead to phenomenological prob¬ 
lems since couplings to the first generation are strongly con¬ 
strained. Therefore, we will assume that the tree-level £ cou¬ 
pling vanishes. This is easily justified in our model because 
this term is not induced via renormalization group running if 
it is zero at some high-energy scale (where one may specu¬ 
late about a ultraviolet completion). The reason is our choice 
of the U(l) x charges li42ll . In addition, we must keep the 1- 
loop induced £ coupling, generated in loops including heavy 
vector-like quarks and leptons, under control. We find 

0 glgX , f m Q \ ™ v , 

Moo r K 16 jt 2 g UJ ' 

Therefore, itiq ~ m R would ensure small E couplings, while al¬ 
lowing large gx . Something required by DM constraints (see 

Sec. [WJ. 


A. Solving the b — t v anomalies 

This model solves the b —► s anomalies in a similar fash¬ 
ion as the one in Ref. Il33ll . The Z' couplings to the SM 
fermions are generated after their mixing with the correspond¬ 
ing vector-like quarks and leptons, as shown in Fig. |T] Ne¬ 
glecting m 2 b < iii 2 q and m 2 <C mj, the resulting Z' couplings 
are found to be 

A „ s = 2 g x X b Q X%v 2 AWI = 2gx|A£| 2 v| 

2m 2 + (|X i e | 2 +|^| 2 )v 2 ’ ^ 2m£+|A£| 2 v 2 - 

( 21 ) 


Therefore, the <f> field will be responsible for the spontaneous 
breaking of U(l)x, giving a mass to the Z 1 , m z i = 2gx\’t/, 
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ultraviolet completion of C\, is obtained with the addition of 
the scalar S = (1,1,0, —4), together with a vector-like (Dirac) 
fermion F = (1,1,0,2). This allows for the Yukawa couplings 
XsSF c F and yLHF, which lead to 0 v after integrating out F 
and L. Although S must get a non-zero VEV in order to break 
lepton number and generate neutrino masses, we note that our 
choice q$ = —4 guarantees that the remnant Z 2 symmetry that 
stabilizes the DM particle % is preserved. 


FIG. 1. Generation of Og and O\o in our model. 


H 


H 


\ / 



1 
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FIG. 2. Non-trivial neutrino mass generation in our setup. 5 is a new 
scalar field with U(l)x charge —4, necessary in order to make the 
operator gauge invariant. 


In case X'q <C 1, the Z! coupling to a pair of SM quarks can 

yb ys* 

be further approximated to A*' ~ gx Q 2 Q ■ These results 

can be combined with Eq. ([8]) in order to determine the al¬ 
lowed ranges for the model parameters that explain the b —> s 
anomalies found by LHCb. 


B. Neutrino masses 

We can further extend the model to get non-zero masses 
for the SM neutrinos. This can be done trivially by adding 
new particles, singlets under t/(l)x, which mediate the stan¬ 
dard mechanisms. For example, the addition of right-handed 
neutrino singlets, Vr = (1,1,0,0), allows for the usual type-I 
seesaw mechanism. 

It is, however, more interesting to consider a mechanism 
that involves the C/(l)x sector of the model. This can be done 
by means of the effective operator 

0 V = \UHH<WS, (22) 

Av 

as shown in Fig. [2] Here S is a new scalar field with with 
U (1 )x charge qs - —4, necessary in order to make the op¬ 
erator gauge invariant. An example model that can serve as 


At high temperatures % attains thermal equilibrium in the 
heat bath through its reactions with the different degrees of 
freedom to which it couples. At decoupling, the dominant 
processes determining the DM yield are the following 2f>2 
reactions 0 

• FI FI' -n- XX* (Higgs portal): processes enabled by the 
hit |H| 2 |5C| 2 coupling. 

• ff f> XX*’ with F standing for the SM and the new 
vector-like quarks and leptons (Z' portal): the former 
enabled by F$m~F mixing and U {1 )x coupling, while 
the latter solely by gx- 

For very heavy DM with m x > m z /, also Z'Z' o XX* can 
take place. Depending on the relative size of these processes 
(pure scalar and Z'-mediated) one can then distinguish several 
scenarios. Among them those solely involving Z'-mediated 
processes, are—arguably—expected to be dominant (they are 
driven by a gauge coupling). Interestingly enough, in that case 
a clear correlation with flavor physics must exist. Note that 
these processes match those in Fig. |T] if one trades one of the 
fermion pairs for XX* ■ Therefore, under the fairly reasonable 
assumption that the Z'-mediated processes play a dominant 
role, an interplay between flavor and DM physics is possible 
establishing. Thus, in turn, further constraining the flavor- 
transition parameters through the restrictions imposed by the 
condition of generating the correct DM relic density (see the 
useful Ref. PTTll f] 

Although our results rely on MicrOmegas l46l . a simple 
analytical discussion is worth doing. The cross section for the 
FF o XX* processes can be estimated to be 


a ( s )' 




8 x 




-f{x f ,x s 


(23) 


where the x, = mj/s and f(xf,x s ) is a kinematic function. De¬ 
pending on the relative size of niy and iuq.l (r = niy/mQx), 
one can distinguish two regimes. For r < 1, DM annihi¬ 
lation processes involve dominantly SM quarks and leptons 
(b/ .si and/.//J. In that case, however, the corresponding cross 


- Note that since m v i > m f . Z' / yj/j plays a subdominant role. 

3 Further constraints such as direct/indirect detection and collider searches 
might be also relevant, see ram- 
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sections are suppressed by chiral/vector-like mixing (see eq. 
(|2T|l), annihilation is rather inefficient and therefore leads to 
an overpopulation of % scalars. For r > 1, processes involving 
the vector-like fields are also possible, and since they are not 
suppressed by mixing parameters they can lead to the appro¬ 
priate DM relic abundance. Note that if the DM and vector¬ 
like fermion mass splitting is small, then the cross section will 
be phase space suppressed: / —> 0. In this case, resonant anni¬ 
hilation is needed to efficiently deplete the scalar % population 
(see Fig. [3J. For large mass splittings, in contrast, annihilation 
is very efficient, and so the correct relic density can be read¬ 
ily obtained. However, reconciling this “scenario” with flavor 
constraints turns out to be tricky. Avoiding the resonance re¬ 
quires my > /n z , which means my m/ q, and a large gx 
close to 1. In that limit. 


✓-.NP 

C 9 


^b 'is 

KqKq 

2 ’ 
m Q 


(24) 


which for typical vector-like quark masses, ihq ~ TeV, im¬ 
plies large chiral/vector-like quark mixing (order one Xq S cou¬ 
plings) in order to be compatible with the LHCb observations. 
This, however, is forbidden by quark flavor constraints (see 
next section). 

In summary, the correct DM relic density can be easily 
produced within our setup provided r ~ 0(1) and DM an¬ 
nihilation proceeds resonantly. Resonant annihilation can be 
avoided for r ^ 1, but finding spots in parameter space consis¬ 
tent with both, DM and flavor constraints seems challenging. 
The right relic density might as well be produced when r < 1, 
but probably this would require being sharply at the Z' reso¬ 
nance. 


V. CONSTRAINTS 


There are constraints on the mass of additional gauge 
bosons coupled to SM states. The ballpark of these limits is 
about 2.5-3.0 TeV for 1/(1) extensions which predict a cou¬ 
pling of the Z' to light quarks of 0(1) l47l . However, in 
our model these couplings are suppressed by the small mix¬ 
ing between the SM quarks and the vector-like states. Thus, 
we nearly get any constraint on the mass of the Z' from LHC 
searches. In addition, we are going to assume in the lepton 
sector only a sizable mixing for muons with the new states. 
Therefore, we get also hardly any limit on my from LEP 
searches. 


A. Collider constraints 

The masses of the vector-like quarks have strong bounds 
from the LHC. Being colored particles, vector-like quarks 
can be efficiently produced in pp collisions, which typically 
pushes their masses towards the TeV scale, see |481 and refer¬ 
ences therein. Since our setup works with vector-like quarks 
with masses at or above the TeV scale, these bounds are easily 
satisfied. 

Vector-like leptons can be searched for at the LHC in the 
standard multi-lepton channels. Using the CMS analysis in 
in, based on searches for final states including 3 or more 
leptons with an integrated luminosity of about 19.5 tb 1 at 
Vs = 8 TeV, Ref. ED obtained a lower limit on the mass 
of doublet vector-like leptons of about 460 GeV, in case they 
decay to electrons or muons, and about 280 GeV, in case they 
decay to tau leptons. More recently. Ref. iBTTl considered the 
analogous ATLAS multi-lepton search l52l . with 20.3 fb 1 
at y/s = 8 TeV, also looking for final states with 3 or more 
leptons. The results were similar to the ones derived from 
the CMS analysis, with lower limits on the mass of doublet 
vector-like leptons of about ~ 500 GeV. 

Finally, one can also derive limits from Z —»• 4£ searches at 
the LHC. However, according to ||33l . these are rather mild, 
only relevant for Z' masses below 100 GeV. 


B. Quark flavor constraints 

Contrary to the model in lf33l . where the muonic cur¬ 
rent contributing to the b —> s observables is purely vectorial 
(£y a £), in our model the current is left-handed and thus con¬ 
tains an axial vector contribution as well (LyaYsV- This is rel¬ 
evant for the B s -A p p decay, especially sensitive to axial- 
vector currents. Currently, there is a little tension between 
the SM prediction for the average time-integrated branching 
ratio, BR(ZL -A p p ), and the CMS and LHCb measure¬ 
ments. While the SM prediction is found to be BR(B S —> 
p + p~) sm = (3.65 ±0.23) x 10~ 9 lf53l , the combination of the 
CMS and LHCb measurements leads to BR(B, -A p l p )exp = 
(2.9 ±0.7) x 10~ 9 (54). In view of this little deficit, the new 
Z' mediated contribution might be potentially welcome. 

The B s —> p p decay is known for its reach in probing 
new physics scenarios (see for example the recent papers lf55l - 
[593 ). In the SM, the amplitude is dominated by the axial 
vector leptonic current, and thus by the C ! ‘ l() Wilson coeffi¬ 
cient (C)q SM ). Similarly, in the model under consideration, 
the main NP contribution is given by a modification of C‘“ () . 
We note that vector contributions vanish for lepton flavor con¬ 
serving channels. Therefore, our prediction for the average 
time-integrated branching ratio is simply given by 


Thus, we can safely assume in the following that the Z' cou¬ 
plings to the first generation of SM fermions can be neglected 
and evade all LEP and LHC limits. Similarly, 2 ! contribu¬ 
tions to flavor observables related to the first generation can 
be safely neglected. Let us now review several relevant con¬ 
straints in our model. 


BR(B S -A p + p ) 


✓tU.SM NP 

L 10 + L 10 

wi,SM 

L 10 


2 

BR(B 5 — >p + p )sm • 
(25) 


Nevertheless, the large experimental error in BR(/i, -a p 1 p ) 
precludes from obtaining useful bounds. Using the SM 
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and the experimental average time-integrated branching ratios 
given above, one finds —0.25 < C^ 0 NP /C^q SM < 0.03 (at the 
la level). By combining C P q NP = —C 9 ' NP with the SM val¬ 
ues for Cg and C p 0 , one finds that in our setup C|q NP /C^ 0 SM ~ 
0.06, which is slightly above the la limit, and thus perfectly 
compatible with the experimental measurement of BR(Z? S —t 
pp ) at the 2 a level. 

Another relevant flavor constraint comes from B s — B s mix¬ 
ing, induced at tree-level by Z' exchange as soon as a non¬ 
zero coupling is considered. Allowing for a 10% de¬ 
viation from the SM expectation in the mixing amplitude, 

|M 12 /M\f - 11 < 0.1, one finds (TT) 

> 244 TeV. (26) 

K s \ 

C. Lepton flavor constraints 

In principle, the mixings with the vector-like leptons can 
induce lepton flavor violating (LFV) processes of the type 
t a —> mediated by the Z' boson. However, one can 

suppress them with a proper parameter choice. Since all b —> s 
anomalies can be simultaneously explained with (only) new 
contributions to Og, X^ ^ 0 is required, but one can choose 
X‘ r ~ = 0. This would eliminate the coupling of the Z' boson 
to electrons and tau leptons, and thus all LFV processes medi¬ 
ated by the Z' boson. This includes LFV in B decays, recently 
suggested in (62) and further studied in ||3 1 1 1631465) . 

Besides the anomalies in b —> s transitions, the LHC might 
have found additional hints for new physics. Recently, the 
CMS collaboration found a 2.4a excess in the h —> xp chan¬ 
nel which translates into BR (h —> xp) = (0. 84 + 037 )% ( 66 ). 
For this result, the collaboration made use of the 2012 
dataset taken at y/s = 8 TeV with an integrated luminosity of 
19.7 fb~ 1 . This large Higgs LFV branching ratio cannot be ac¬ 
commodated in our setup, since vector-like leptons are known 
to be unable to reach such LFV rates in Higgs decays due 
to limits from the radiative x —> py decay EQ). However, as 
suggested in I67ll68l and recently confirmed in several works 
(Ml HIM], a simple extension with a second Higgs doublet 
suffices to explain the CMS hint. 

Finally, we comment on limits from the (former) non¬ 
observation of lepton universality volation. Indeed, very 
strong bounds have been derived in many scenarios due to 
the non-observation of lepton universality violating effects in, 
for example, pion and kaon decays (see for example (724)7410 . 
However, these are absent in our model due to the suppression 
of the vector-like and first generation SM quark mixing. 

D. Precision measurements 

There are several precision measurements in the lepton sec¬ 
tor which might be potentially sensitive to the Z! interactions 
considered in this model. 

First of all, there is the so-called neutrino trident produc¬ 
tion (60) . This is the production of a muon anti-muon pair by 


scattering of muon neutrinos in the Coulomb field of a target 
nucleus. The cross-section gets additional contributions in the 
presence of a new Z' boson that couples to the muons, and 
these can be used to constrain A+“ and niyj. The impact of this 
bound on Z' models designed to solve the b —> ,v transitions 
was studied in 13311341 . Using the CCRF measurement of the 
trident cross-section ED, one can derive bounds on the mass 
of the Z' boson and its coupling to muons. For scenarios with 
left-handed Z' coupling to muons, one finds m 

TXWi ~ 470GeV. (27) 

K\ 

Finally, there are two other constraints related to precision 
measurements in leptonic observables: (g — 2)„ and the Z bo¬ 
son couplings to the muon (which are modified at the 1 -loop 
level due to the Z' interactions). However, in (33) these are 
found to be relatively mild. For this reason we will not con¬ 
sider them in our numerical analysis. 

VI. NUMERICAL RESULTS 

For the numerical analysis we have implemented the model 
into the Mathematica package SARAH (7614801 and generated 
Fortran code to link this model to SPheno (81)1821 . The big 
advantage of this setup is that we can cross-check our tree- 
level approximations with a full numerical calculation. More¬ 
over, based on the FlavorKit interface [83), the generated 
SPheno version does not only calculate the tree-level contribu¬ 
tions to the Wilson coefficients but also all 1-loop corrections. 
This way, we can see in what parameter range the tree-level 
approximation works well and where it breaks down. In gen¬ 
eral, one can expect loop effects to become important when 
the mixing between SM quarks and the vector-like states is 
sizable: this leads to significant changes in the SM-like box 
and penguin contributions to Cg and C io- This alone would 
not be necessarily a problem. However, the same mixing in¬ 
fluences also loop contributions to the other Wilson coeffi¬ 
cients, in particular C 7 and Cg. These effects have to be small 
in order to satisfy the results of global fits to b —> s data. 

As second step, we used the CalcHep (841 1 85) output of SARAH 
to implement the model in MicrOmegas to calculate the relic 
density. Interfacing the parameter values between SPheno and 
MicrOmegas is done via the SLHA+ functionality of CalcHep 
& 

Our main numerical results are summarized in Figs. [3] and 
[4] In Fig. [3] we show the dark matter relic density in the 
(gx , niyj ) plane, together with the ratio C^ P /C| M . The other 
important parameters entering this calculation have been fixed 
to 

Xq = X S Q = 0.025, XI = 0.5 
niQ = mi = 1 TeV, m \ = 1 TeV 2 

First of all, we see that with these parameters the model is 
perfectly compatible with the constraints discussed in Sec¬ 
tion [V] Furthermore, there is a region for moderately large 
gx — 0.3 where the DM constraint can be satisfied and Cg = 
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FIG. 3. Contours for constant C^ P /C| M and log(£^DM^ 2 ) (dashed 
black) in the (gx, m z') plane. For C^ p /Cg M we show the full 1-loop 
results via the black lines, while the dotted grey lines give the values 
using the tree-level approximation only. 



FIG. 4. Contours for constant Cg^/Cg 1 ^ (full black lines) and 
Cf/Cf 1 (dotted grey lines) in the (gx,mz') plane. 


— 12% x Cg M holds. To reduce the relic density to £llr ~ 0.1 
one has to be rather close the the Z' resonance around 2 TeV. In 
the same Figure we show also Cg p when using the tree-level 
approximation only. In the interesting region where the fla¬ 
vor anomalies can be explained, this approximation is work¬ 
ing quite well. However, for decreasing gx one would expect 
smaller values for C^ p at tree-level, while we find a large in¬ 
crement. The reason is that for constant m z i and smaller gx, 
v,|, is increasing, leading to a larger mixing in the quark sector 


with the vector-like states. Thus, loop effects become impor¬ 
tant and quickly dominate the behavior. These effects are not 
only expected to show up in C^ p but also in the other Wilson 
coefficients as stated above. To demonstrate this, we plot in 
Fig. [rj the contours for Cq 9 /C$ M and C^ p /Cf M in the same 
plane. Obviously, we find a similar enhancement for small 
gx for both coefficients. However, for the most interesting 
region where tree-level contributions to C^ p dominates and 
can explain the anomalies, the change in C 7 is very moderate: 
C 7 np /C 7 sm < 1%. 


VII. SUMMARY AND DISCUSSION 

We have presented a DM model that successfully addresses 
the anomalies in b —► s transitions recently found by the LHCb 
collaboration. In our setup, B mesons decay through the Z' 
into the dark sector, which mixes with the SM via Yukawa 
couplings. An interesting connection between DM and flavor 
physics emerges, leading to a very constrained scenario where 
the NP masses and parameters are restricted to lie in thin re¬ 
gions of the full parameter space. 

The model offers several interesting phenomenological 
possibilities. Due to the mixing between the SM muons and 
the vector-like lepton doublet, our model predicts a reduction 
of the h — /j — /j coupling, ghuu, with respect to the SM predic¬ 
tion. The most direct probe of this coupling is, of course, the 
Higgs boson decay into a pair of muons, h -4 /j 1 /r , experi¬ 
mentally very challenging in a hadron collider due to the ex¬ 
pected low branching ratio (0(10 4 ) in the SM). In 2013, the 
ATLAS collaboration presented results based on an integrated 
luminosity of 20.7 fb -1 in collisions at t/s = 8 TeV, without 
any evidence of a signal ED- For in/, = 125 GeV, they ob¬ 
tained an upper limit on the qq signal strength, /./ H „ <9.8 (at 
95% CL). Assuming a SM-like Higgs boson production, this 
limit translates into gi 1/Jf , < 3.1g™. Obviously, our model is 
well within the limit, since it predicts a reduction of the cou¬ 
pling. Nevertheless, it is worth pointing out that a precise de¬ 
termination of this coupling in a future linear collider would 
be a strong test of the model |^] 

Let us also comment on the mixings between the vector¬ 
like fermions and the SM ones. These are determined by free 
Yukawa parameters which, in principle, are naturally expected 
to be of the same order for all SM generations. However, in 
order to avoid conflict with existing experimental data, one is 
forced to strongly suppress the mixings with the first genera¬ 
tion, implying X d Q,X e L -C X S q ,X^ x . This might seems like an 
ad-hoc assumption. However, it just reflects the constant need 
for a theory of flavor, also required to understand the Yukawa 
structure of the SM0 


4 According to HD, CLIC would be able to measure the h — /j — n coupling 
with an uncertainty of about 19%. 

5 We note that although the lepton sector of (33) naturally explains the re¬ 
quired hierarchy for the lepton couplings, in the quark sector the situation 
is exactly the same as in our model. 
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Finally, one can envisage several phenomenological direc¬ 
tions in which this model can be further explored. By allowing 
for non-zero couplings to first generation quarks and leptons, 
the impact of the model gets extended to many additional ob¬ 
servables of interest. In particular, the LHC would be able 
to probe the parameter space of the model via the produc¬ 
tion of the Z' boson. Furthermore, since the scalar DM car¬ 
ries, at early times, a conserved U(l)x charge it will develop 
an asymmetry, provided chemical equilibrium with the ther¬ 
mal bath is guaranteed. In such a case (asymmetric DM sce¬ 
nario f88l ), a link between the baryon and DM asymmetries 
is something worth exploring. 

NOTE ADDED 

An update of the global fit fTTI was recently presented in 
l89l . While this would change slightly the numerical output 


of our analysis, our conclusions would remain unchanged. 
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